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The strong time-reversal symmetric (TRS) topological insulator (TI) in three space dimensions 
features gapless surface states in the form of massless Dirac fermions. We study these surface 
states with the method of bosonization, and find that the resulting bosonic theory has a topological 
contribution due to the parity anomaly of the surface Dirac fermions. We argue that the presence of 
a quantum anomaly is, in fact, the main reason for the existence of a surface state, by the principle 
that anomalies of a surface and bulk must cancel. Inspecting other classes of topological insulators, 
we argue that this principle holds in general. Moving beyond purely topological considerations, we 
incorporate the dynamics of the surface electron states into the bosonic theory. Additionally, we 
discuss the thermodynamics of the bosonic theory and propose a representation of the surface Dirac 
fermions in terms of the bosonic fields. 
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I. INTRODUCTION 

Topological insulators [IH3] are a recently discovered 
class of materials whose most salient feature is that they 
are insulating in the bulk but feature gapless, topolog- 
ically protected surface states. Here, we will focus on 
the so-called "time-reversal symmetric (TRS) topologi- 
cal insulator (TI) in three space dimensions" 0t3 ; which 
is protected against time reversal invariant perturba- 
tions. In the regime of low energies, the surface states 
of the TRS topological insulator in three space dimen- 
sions are described by 2 + 1 dimensional massless Dirac 
fermions [§], Q , which obey the massless Dirac equation 
7 M 9 M V' = 0, fJ, = 0, 1, 2. Within this fcrmionic approach, 
the effect of electron interactions on the magnetic sus- 
ceptibility and on magnetic symmetry breaking has been 
studied [lMl- 

In this article, we discuss a description of TI surface 
states in terms of a bosonic field theory, which contains 
both a universal topological term and a non-universal 
term describing the electron-hole continuum of massless 
Dirac electrons. The idea of describing fcrmionic systems 
in terms of composite bosonic excitations — like the elec- 
tron density — is known as bosonization and has been 
successfully applied to 1 + 1 dimensional systems fl4r - 
[l6| ; here we argue that it is useful in 2 + 1 dimensional 
systems as well. 

More specifically, we will use functional bosonization 
[lll-IHl to construct a topological field theory [H HI 
that connects bulk and surface states in a satisfying man- 
ner. We emphasize the important role of the parity 
anomaly 23l425| and of quantum anomalies in general 
l3~l| for understanding the existence and topological 
nature of surface states, and show how the dynamics of 
the surface states can be incorporated into the field the- 
ory. By construction, the resulting bosonic field theory 
correctly describes the electromagnetic response in the 
long wavelength limit, but reproduces the fcrmionic ther- 
modynamics only qualitatively. We explain that func- 



tional bosonization cannot be expected to quantitatively 
reproduce thermodynamics quantities from the bosonic 
effective action alone, and that the difference between 
fcrmionic and bosonic thermodynamics is accounted for 
by normalization factors of auxiliary functional integrals. 
Moreover, we discuss the current status of the refermion- 
ization^^, l32l l33j program in 2 + 1 dimensions, which 
aims at studying the underlying fermion degrees of free- 
dom in terms of the bosonic theory and whose main at- 
traction is that it could simplify the study of strong elec- 
tron interactions. 

The structure of the paper is as follows: First, in Sec- 
tion nn we discuss the quantum Hall effect as an exam- 
ple for demonstrating how the existence of surface states 
follows from the principle that quantum anomalies of 
both surface and bulk must cancel, in this case the chiral 
anomaly. Then, we apply this idea to the TRS topolog- 
ical insulator in three space dimensions and argue that 
here the parity anomaly is responsible for the existence 
of surface states. In Section IIII1 we review the method 
of functional bosonization and apply it to the TRS topo- 
logical insulator to obtain a topological field theory that 
connects bulk and surface via the parity anomaly. Af- 
ter that, in Section IIV1 we focus on the surface states 
and incorporate the electron dynamics into the bosonic 
field theory by considering the non-topological parts of 
the effective electromagnetic response. In Section [Vj we 
calculate the partition function of both the bosonic and 
fcrmionic surface theory, and discuss the discrepancy in 
the specific heat between the two approaches. Finally, in 
Section |VT1 we review and use the so-called tomographic 
transform to propose a representation of the electron field 
in terms of the bosonic theory. A summary I VIII concludes 
the main text. Appendix lAl collects some tools for ma- 
nipulating the surface Chern-Simons term. Finally, in 
Appendix [B] we return to quantum anomalies and dis- 
cuss how the quantum spin Hall effect relates to the axial 
anomaly. 
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II. CANCELLATION OF QUANTUM 
ANOMALIES 



Our goal is to describe topological insulators by a field 
theory of bosons that arise as collective excitations of 
the underlying electron field. This approach has been 
very successful for studying the integer and fractional 
quantum Hall effects [34| . 

Such a field theory is necessarily topological in nature 
and its defining feature must be that it correctly pre- 
dicts the topologically protected surface states. We will 
now argue that the main mechanism for the emergence 
of surface states is the existence of a quantum anomaly 
of the bulk states together with the requirement that 
this anomaly is cancelled by the surface states. In par- 
ticular, for the TRS topological insulator in three space 
dimensions, the quantum anomaly in question is the par- 
ity anomaly; for other topological phases, the relevant 
anomalies are different. 

In quantum field theory, the notion of a quantum 
anomaly (35j refers to the phenomenon that a symme- 
try, like gauge invariance or time-reversal symmetry, may 
be present in the classical Lagrangian, but is lost once 
the theory is quantized. In particular, this means that 
the response of the quantum system to an external field, 
like the electromagnetic field, docs not have the expected 
symmetry cither. We will now describe the response of 
some classes of topological insulators and see that the in- 
troduction of a boundary breaks a symmetry of the bulk, 
which should be restored by the quantum anomaly of an 
emerging surface state. 

As a warm-up, let us first consider the prototype ex- 
ample of a topological phase, namely the two-dimensional 
quantum Hall system. Given an external electromagnetic 
gauge field A^, the real time effective action of the bulk 
is of Chern-Simons type (3~H 

S eS [^] = ~- J d 2 xdt e^ x A^A x (I) 

and describes a Hall response with Hall conductivity a xy . 
For periodic boundary conditions, this action is invariant 
under a gauge transformation A^ — > A^ + <9 M A, which 



implies that the current = is conserved, d^j^ 

I i as expected. 

However, if we restrict this effective action to a domain 
with boundary, for instance to the half-space x\ > 0, we 
see that the corresponding action 

Soft, buik[AJ = J d 2 xdt 6(x 1 )s^ x A fi d 1/ A x (2) 

gives a current that is no longer conserved 



Wbulk — 



SA L 



= -^S( Xl )e l ^d v A x . (3) 



This is not unexpected: a Hall current flows towards the 
edge of the sample, but cannot continue in the vacuum, 



so there is a defect at the edge. Of course, it is a general 
physical principle that the total current must be con- 
served and we conclude that there must exist an edge 
state that compensates the loss [HI, 



M _ 



^Jcdgc 



-fS( Xl )e 



d v A> 



+ ie A dge) = 0- 



(4) 



(5) 



Note the peculiar fact that the fermionic theory at the 
edge individually also violates current conservation. This 
is the chiral anomaly [3l| in I + I dimensions, and the 
compensation presented here is the Callan-Harvey mech- 
anism [37j • Of course, the response of the total system is 
conserved, it is just the individual parts that are anoma- 
lous separately. This means that the anomalous part of 
the response of the edge state is not observable] if any- 
thing, its only observable consequence is the existence of 
the edge state in the first place, and of course any normal 
responses that follow from that. 

We now present a similar argument for the surface 
state of the strong TRS topological insulator in three 
space dimensions. This time, however, the surface states 
are governed by the parity anomaly [23l [24j and the phys- 
ical symmetry principle we invoke is not current conser- 
vation (which follows from gauge symmetry), but time- 
reversal symmetry. 

The bulk electromagnetic response of a TRS topologi- 
cal insulator in three dimensions with periodic boundary 
conditions has been derived KlM ulSl as 



e 2 

S oS [A»] = I d 3 xdt E • B 

27T he 



8tt 



J d 3 xdt s^^d^dpAv 



(6) 
(7) 



where we have absorbed [59( the electron charge e into the 
definition of the gauge field A^. Displayed is the topolog- 
ical part of the response, the dots indicate that there may 
be additional, non-topological contributions. Parameter 
values 9 = ±7r indicate a topologically non-trivial insu- 
lator, E and B are the electric and the magnetic field 
respectively. It is not obvious why this action should be 
time-reversal invariant, since after all, time reversal maps 
E • B — > — E • B. However, for periodic boundary con- 
ditions, the integral / d 3 xdt e^ pa dpA v d p A a is actually 
quantized to an integer multiple of 87r 2 , this is a general 
property of gauge connections related to the so-called sec- 
ond Chern class Jijjj]. Thus, time-reversal does change the 
action by AS = 28n = 27rn, but only the Feynman path 
amplitude e has physical meaning in quantum theory, 
and this one remains invariant. 

Once again, let us now restrict the action to a domain 
Q with boundary dQ. The bulk action can be written as 
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a boundary integral 



SbuikK] = A / d3xdt e^d^dpA* (8) 



8tt 2 



<«2 



(9) 



which is no longer quantized or time-reversal invariant. 
But since the total action of the system must be time- 
reversal invariant, there must exist a surface state whose 
response makes the total action S^uik + Surface time- 
reversal invariant, for instance by 

Surface = - A / d2xdt £ Upa A v d p A a + .... (10) 

8tt 2 J an 

This effective action precisely matches the parity 
anomaly [HI, [24| of the 2+1 dimensional massless Dirac 
fcrmion, which is the surface state that we would ob- 
tain by solving a fermionic (model) Hamiltonian directly. 
Again, this is the topological part of the action, the dots 
indicate that other, time-reversal invariant contributions 
are possible, for instance those that describe the dynam- 
ics of the surface state. We will discuss them in Sec- 
tion ED 

This topological contribution of the surface looks like a 
half-integer quantum Hall effect on the boundary. How- 
ever, we stress that it cancels out in the total electromag- 
netic response and cannot be observed experimentally. 
This was recently also shown by an explicit computa- 
tion of the fermion Green's function in the presence of a 
boundary [23]. In other words, the bulk action (|6|) cor- 
rectly describes a TRS topological insulator with gapless 
surface states, in combination with the anomaly cancel- 
lation mechanism presented here. There is no need to 
artificially break time-reversal symmetry in order to de- 
scribe TRS topological insulators. 

In general, we believe that this mechanism of anomaly 
cancellation applies to all topological insulators with any 
symmetry and in any dimension. The gapless states on 
the boundary should always feature a quantum anomaly 
that cancels with the bulk; this is why they exist in the 
first place. Recent theoretical evidence |25l - l31| has estab- 
lished a general connection between topological insula- 
tors and quantum anomalies. Importantly, all symmetry 
classes 0] of topological insulators are related to anomaly 
polynomials [26[. Our conjecture about the connection 
between boundary states and quantum anomalies is sup- 
ported by the following collection of examples: 

i) Integer quantum Hall effect in 2 space dimensions, 
symmetry class A. Transfer of charge between bulk and 
edge due to the chiral anomaly of the boundary states. 

ii) TRS topological insulator in 3 space dimensions, 
symmetry class AIL Cancellation of the parity anomaly 
on the boundary as discussed above. 

hi) Quantum spin Hall insulator in 2 space dimensions, 
symmetry class AIL Combined spin or charge flow be- 
tween edge and bulk due to the axial anomaly [4l| . We 



delegate the discussion of this instructive example to Ap- 
pendix [B] as we want to focus on three space dimensions 
in the main text. 

iv) Particle-hole symmetric topological superconductor 
in 2 space dimensions, symmetry class D. The boundary 
states are described by Majorana fermions, but these do 
not carry charge, so they cannot give rise to a quantum 
anomaly involving the electromagnetic field. The bound- 
ary modes do, however, conduct heat and it has recently 
been shown j27| that the (thermal) energy and momen- 
tum current has a quantum anomaly and flows from the 
bulk into the edge. As energy and momentum are the 
conserved quantities corresponding to spacetimc transla- 
tions, they can be probed by distorting spacetime, i.e. by 
introducing a gravitational field [iiOj Hence, it is a gravi- 
tational anomaly [27], Hi], S3 that explains the emergence 
of Majorana fermions at the boundary. 

We believe that this list of examples is sufficiently di- 
verse to support our conjecture about the connection be- 
tween quantum anomalies and topological surface states, 
as it includes the basic known quantum anomalies, fea- 
tures different cancellation mechanisms and mentions 
both primary insulators and secondary descendants ob- 
tained by dimensional reduction [38j . 



III. TOPOLOGICAL FIELD THEORY 

Having understood the parity anomaly as the topo- 
logical mechanism that requires the existence of surface 
states, we can now derive topological field theories for 
surface and bulk. This can be done with a general 
method called functional bosonization [T341H > which we 
will now review to make our discussion self-contained. 

The main idea of functional bosonization is to start 
with some effective action S' c ff[A A1 ] for the external vec- 
tor potential that was obtained by integrating out the 
fermion degrees of freedom. Then, we subject this action 
to a so-called duality transform [17J , which yields bosonic 
fields that generate the same effective action. These fields 
can then be interpreted as collective charge excitations 
of fermions. In a sense, we "undo" the process of inte- 
grating out the matter fields. 

For concreteness, let us imagine the example of a mass- 
less 2 + 1-dimensional Dirac fermion, which represents 
the surface state. Assume for a moment that we have 
somehow managed to calculate its partition function and 
effective action 



Z[A, 



DipDi/jexp i / d 2 xdt V^7 M (d M - iA^ip 



when coupled to an external gauge field A^. We know 
that this functional is gauge invariant, so we could as well 
average over all configurations that are connected by a 
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gauge transformation 



Z[A^\ = J DA Z[A^ + dpA] 



(12) 



up to a normalization factor that we ignore for now. We 
can write this as an integration over a field a p that is con- 
strained to be a pure gauge a p = d^A. By the Poincare 
lemma, this constraint is equivalent to the vanishing of 
the curl e p ' J> 'd v a\, which we can impose with a delta 
function 



Da. 6(e^ x d„a x )Z[A p + a p 



(13) 



Shifting a p — > a p — A p and enforcing the delta function 
constraint in terms of a Lagrange multiplier Wn , we ob- 
tain 

Z[A„] = J Da^DW^ exp (iS eS [a„] (14) 

" h / Sxdt e ^ Xw ^ a ^ - 

=: J Da^DWf, exp (iS hoson [a^, W M , A^]j . (15) 

The expression in the exponential is the desired action for 
the bosonic fields a p and W^. The whole procedure is es- 
sentially just a fancy way to rewrite our starting point, 
the effective fermionic action, in terms of the charge cur- 
rent 

= SS h [a,W,A] s 1 ^ 

27T 

which has been expressed as the curl of a field because 
it needs to fulfill = 0. 

Applying this method to the analysis of surface states, 
i.e. to the effective action ([TO)) corresponding to the parity 
anomaly, we obtain a bosonized theory 



Air 



- £^W^a p + e^"W^A p \ . (17) 

This is the topological field theory of the surface state; it 
is equivalent to the celebrated Chern-Simons theory [3J] . 
It is topological as the corresponding equations of motion 
for the bosonic fields are devoid of dynamics. Of course, 
the actual massless Dirac fermions on the surface do have 
dynamics, i.e. their effective actions contains terms be- 
sides the parity anomaly. We will discuss these shortly, 
but for now let us focus on the topological aspects fol- 
lowing from the parity anomaly. 

In the same way, a topological field theory for the 
bulk can be obtained j20j . In three space dimensions, 
the derivation just presented has to be changed slightly: 
the Lagrange multiplier needs to incorporate more con- 
straints to make the field a M a pure gauge, and so will 



become a tensor field & M „ instead of a vector field W^. 
Starting with the effective bulk action <j6j> , one obtains 
the topological field theory [20l | 



SbvMt[au,buv,Au] = — / d xdt 
27T Jn 



47T 



e^ pa b^d pa<J + e^ pa b pv d p A a 



(18) 



for the bulk. 

We can now connect the topological field theories for 
bulk and surface. Cancellation of the parity anomaly 
dictates that the topological bulk and surface actions 
should cancel, so we have to be able to express the former 
as an integral over the boundary. This is indeed possi- 
ble: First, the term e pvpa d^AydpCLa can be moved to the 
boundary similar to Equation Second, we can utilize 
the Euler-Lagrange equations from variation of the field 
a p to obtain the condition e^^d^byx = 0. This means 
that the Lagrange multiplier is forced to be a pure gauge 
bpu = —^{dpWv — dvWp). Using this relation in Eq. (jTSf . 
we can move the remaining two terms to the boundary as 
well and we obtain precisely the negative of the surface 
action. 

Compared to previous proposals for field theories of 
topological insulators [2(| HH , our approach is guided by 
the principle of anomaly cancellation between bulk and 
surface states. In contrast to [22[, we arrive at a surface 
action which is not invariant under time reversal on its 
own, but together with the bulk action is time reversal 
invariant, and elucidates the topological meaning of the 
parity anomaly. 



IV. SURFACE DYNAMICS 

While the cancellation of quantum anomalies captures 
the topological content of the field theories for surface 
and bulk, it is clear that topology alone cannot give any 
information about the dynamics of the surface states. 
For instance, the edge action for the quantum Hall ef- 
fect may contain any additional term that is compatible 
with gauge invariance [3(|. Likewise, the TRS topologi- 
cal insulator in three space dimensions may contain terms 
compatible with both time-reversal invariance and gauge 
invariance. Unlike previously suggested [22| , these terms 
are not restricted to represent a potential energy. In fact, 
we will obtain kinetic terms that contain both space and 
time derivatives, even in a non-local fashion. 

To obtain the dynamic parts of the surface theory, we 
have to start with the full effective action for the 2 + 1 
dimensional massless Dirac fermions on the surface. So 
far, this effective action could not be evaluated exactly, 
but was only computed to quadratic order in the field a M 
[43M45| ] . Using this result in the functional bosonization 



5 



formula ([Ml), we obtain the full surface action 



of Euclidean momentum q as 



S B w* BC e[a li ,Wit,A l t] = J d 2 xdt 



8tt 2 



+ l^(- n) V2 a T + 0(a 4 ) 
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~e^\a^-A^)d v Wx 

Z7T 



(19) 



where denotes the transversal components of the 
vector field a M . The first term is the by now familiar 
parity anomaly, while the second one describes the elec- 
tron dynamics. The dynamic term looks quite unusual, 
though: It is not a purely potential energy, for it involves 
the temporal derivatives of the field. It is not the 
Maxwell term either, because it features the square root 
of the d'Alembert operator (— D) 1 / 2 , a non-local pseudo- 
differential operator. Fortunately, the retarded Green's 
function G R of this operator is well defined in 2 + 1 di- 
mensional space time (46j 



G R (x,x') 



1 



An x — x 



T At-t'-\i 



(20) 



This is just the 3 + 1 dimensional wave propagator re- 
stricted to the plane X3 = 0. Note that only values on 
the boundary of the light cone t — t' = |x — x'| con- 
tribute to the propagator; in mathematics, this property 
is sometimes called Huygcn's principle. Compare this 
to the usual 2 + 1 dimensional wave propagator, which is 
logarithmic and also involves points inside the light cone, 
not just the boundary. 

In momentum space, in particular in Euclidean space- 
time, the square root term is not unusual either, it 
is simply the magnitude of the Euclidean momentum 



(-□) 



1/2 



\q\. This dynamic term corresponds to the 



density-density correlation function [4£ 



(j°(-q)3°(q)) 



(21) 



which yields \q\ when combined with the projection onto 
the transversal components a^(q) = (5^ — ^^-)a ll (q). 

At this point, we believe it is useful to discuss the effec- 
tive action for massive 2 + 1 dimensional Dirac fcrmions 
as well, since we feel that there has been some confusion 
in the literature about the implications of time-reversal 
symmetry breaking via a Zeeman field on the surface 
states. For 2+1 dimensional fermions, the coupling to a 
Zeeman field has the same form as a mass term, and we 
will denote the magnitude of the Zeeman energy by M 
in the following. In Euclidean spacetime, the quadratic 
part of the action has been calculated [43, |44[ in terms 



-^e^xqx + -^e^xqxh 



(K 

\2M 



, (z wA M t f \q\\ 



A„(q) + 0(A 3 ) (22) 



where M is the fermion mass, fi(x) = tan _1 (a;)/|x| and 
f-2{x) behaves like x for small x <C 1 but approaches 1 
for large x S> 1. The first term is the familiar parity 
anomaly, the second term is an additional surface Hall 
effect, and the final term describes the fermion dynam- 
ics. The mass M of the fermions sets an energy scale. 
For momenta small compared to the mass \q\ <C M, the 
dynamics resemble a Maxwell term with a coefficient . 
In addition, a surface Hall effect with a Hall conductivity 



a 



1 M 
X V 2 \M\ h 



is generated, which corresponds to the zero 
energy Landau level of relativistic fermions. However, 
for momenta large compared to the mass, \q\ 3> M, the 
square root of the wave operator (— D) 1 / 2 = \q\ dominates 
the dynamics and the Hall-like term vanishes. Whether 
a Hall effect can be measured in a static charge trans- 
port experiment [1, [4tJ is now a delicate question as the 
static limit \q\ — > must be considered in relation to the 
scales given by the Zeeman magnetic field M and the 
temperature. In any case, we want to emphasize that 
the surface Hall effect induced by the mass is indepen- 
dent of the parity anomaly, as the latter has its origin 
in the need to regularize at large momenta and is not 
observable, because it cancels with the bulk. 

It seems somewhat unsatisfactory that no educated 
guess could have predicted the unusual term describing 
the electron dynamics from the topological field theory 
(fT7]l alone. (In contrast, for the quantum Hall effect, one 
can add a simple potential energy or use a generalized 
gauge fixing to get the right edge dynamics.) However, 
it has been noted [l!| that the parity anomaly subsumes 
the dynamic term, at least in the following sense: Any 
term that is gauge invariant can be absorbed into the 
Chern-Simons term by means of a non-linear field redefi- 
nition. Let us illustrate this for our particular term that 
describes the dynamics. Define a new field a M such that 



~T 

% 



cos/3 ■ al + sin/3 • e^d^-U^^al (23) 

(24) 



with /3 a complex number such that sin(2/3) = Using 
the calculations in Appendix [XJ we see that the Chern- 
Simons term for this new field will expand into the parity 
anomaly plus the dynamic term for the old field 



7. c^^fl n — J!—e^ X n r) n -X- 1 n T ( Fl\ 1 / 2 n T 
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(25) 
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Even beyond second order in perturbation theory, it is 
always possible to find a field a M depending in a non- 
linear way on such that the action becomes the Chern- 
Simons term for the a M field In this sense, the sur- 
face action is given by a single Chern-Simons term arising 
from the parity anomaly, at least modulo a field redefi- 
nition. However, we note that such a redefinition, being 
non-linear in general, gives rise to a non-trivial Jacobian 
in the path integral. 

The bosonic action presented here may also provide 
a good starting point for studying the effect of electron 
interactions. After all, the electron density is given by 
j° = l/27r(9iW / 2 — diW2) and we can model for instance 
the Coulomb interaction by adding a term 

ScouWbtW/J = ^ J dtd 2 xd 2 y (d Xl W 2 - d xa Wi)- 

.—L—fdy^-dyM (26) 

to the bosonic action. Integrating out the field gives 
an action that is again quadratic in the field and 
hence exactly solvable. Of course, this assumes that we 
expand the effective electron action only to quadratic or- 
der, i.e. that we stick to a random phase approximation. 



SURFACE THERMODYNAMICS 



calculation [43, |49[ yields the internal energy 



While the bosonized action ([19]) describes the low- 
energy dynamical and topological properties of the sur- 
face electrons correctly, we would like to remark that — 
unlike in 1 + 1 dimensions — thermodynamic quantities 
like the specific heat are not captured by this action. Af- 
ter all, the field W M is just a Lagrange multiplier and 
the field a M is only a gauge degree of freedom, and hence 
their partition function is trivial from a thermodynamic 
point of view. Of course, the total partition function still 
contains the fermionic thermodynamics, but one needs to 
carefully keep track of prefactors in the bosonization pro- 
cedure to sec this. In contrast, in 1 + 1 dimensions these 
prefactors can be interpreted as a fermionic determinant. 
As a consequence, the fermionic and bosonic action give 
rise to the same free energy. 

For completeness, let us present the calculation in some 
detail. First, we consider a masslcss Dirac fcrmion. The 
logarithm of its partition function can be obtained from 
the functional field integral [HI as 



lnZf c 



1 



lndet( 7 ' i a A1 ) = -lndct(-a i 



Trln(^lc^) 
(27) 

where f3 = 1/kT is the inverse temperature, q — (<7o,q) 
denotes Euclidean momentum and = = 
is a fermionic Matsubara frequency. In the last equality 
in Eq. ([2"T[) . the momentum \q\ 2 has gained a factor of 
P 2 due to the change of the path integral measure under 
Fourier transformation of the field variable. A standard 



d_ 

d/3 

VT 



\aZ = 2V 



3 C(3)3 



d 2 g 

(2tt)2 



1 



71 2 

and hence the specific heat 



CV.fermi 



1 

V 



dU_ 
df 



3T 



2 C(3)3 



(28) 



(29) 



Now, let us compute the partition function for the 
bosonic action ([19p . To this end, we have to express it in 
terms of the Euclidean momentum q, similar to Eq. (|22| . 
Since this will prove useful again later, we have moved 
the corresponding formula (jA7[) to the Appendix [A] Ul- 
timately, the result is that using two matrices P~^ v , P~ v 
that project on the transversal degrees of freedom, we 
can conveniently express the action as 



(30) 



where 



DM = 

±((l + ia)P+ + (l-ia)P-) -±(iP+-iP- 
-£{iP+-iP-) 



(31) 



is a quadratic kernel, and a = ^f- is a constant related to 
the prefactor of the Chern-Simons term. Now, we have to 
project onto the physical transversal degrees of freedom, 
as the matrix would be singular otherwise. Choosing a 
convenient basis, we can write the projections as 2 x 2 
matrix blocks to obtain 



dct D t _ iu {q) = det \q\ 



10 1 



+ ia) 





_ i 

4-7T 


\ 





i(l - ia) 
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-Itt 





0/ 



(4tt) 4 



(32) 



We note that the precise numerical value of the 4x4 
determinant does not matter as long as it is independent 
of temperature, since it contributes only additively to the 
free energy. The important result is that we have four 
bosonic degrees of freedom with a kernel proportional to 
\q\. From this follows the partition function 



\nZ hom = iTrln(/3 4 M 4 ) 
i i 



const. 



Trln(/3 2 | 9 | 2 ) + const. 



(33) 



which corresponds to two copies of the canonical bose 
gas. Once again, the internal energy is calculated by 
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standard methods to be 



t/bc 



d_ 

"dp 



In Z = 2V 



d 2 q 



1 



2 f%l + 1 



yji3 C(3) , 



(34) 



and this result apparently does not match the internal 
energy calculated for the fermions ([28]) . 

It may be useful to compare the above calculation to 
that for the Maxwell-Chern-Simons Theory [Ho||, whose 
action is similar to the action (|19p except that our non- 
local square root (— D) 1 / 2 is replaced by a standard 
Maxwell term l/2aj(— □)aj- The main difference is that 
the Maxwell term is governed by momentum squared q 2 , 
while the Chern-Simons term is proportional to the ab- 
solute value \q\, so that the form of the determinant is 
det(g 2 + K\q\) = det(|g|) det(|q| + K), where the second 
factor can be interpreted as arising from a massive bo- 
son. In our case, however, both the dynamics and the 
Chern-Simons term are governed by |g|, so that the sur- 
face bosons are related to massless canonical bosons. 

However, the apparent discrepancy between the 
fermionic and bosonic specific heat is resolved when tak- 
ing normalization factors of auxiliary functional integrals 
into account. When dealing with the longitudinal degrees 
of freedom, it is important to keep track of Fadeev-Popov 
ghosts, i.e. functional determinants arising from gauge 
fixing [48[ . But since the field a M is only a gauge degree 
of freedom and a Lagrange multiplier, the ghost con- 
tributions will conspire to cancel precisely the partition 
function calculated above. To see this in detail, we have 
to go back to the functional bosonization procedure (fT4|) 
and keep track of prefactors. In particular, averaging 
over the gauge freedom becomes 

Z[A^]=N A J DA Z\A il + dp A] 

= N A J DADa„ S(a„ - d^A)Z[A^ + a») 

= N A J DADa^ S(A - g) det(ft)- 

• 5(d a ai - diao)8{d 2 a - d a 2 )Z[A f _ l + o M ] 



= -^A / det(Sb)- (35) 

• S(d ai ~ dia )S{d 2 a Q - d a 2 )Z[A f _ l + a M ]. 

and expressing the constraints with Lagrange-multipliers 
gives 

Z[A^]=N A N 1 N 2 J Da f ,DW 1 DW 2 dct(d )Z[A^ + a< J- 

d 2 xdt s^ x Wjd v a x ^ 
N A N W I Da^DW^ 5(W )det(d )Z[a fl }- (36) 



■ exp ( " h I d2xdt e ^ Xw n d »( a * - A * 



The constraint on Wo is a gauge fixing condition and 
corresponds to a Fadeev-Popov ghost |48| 



S(W )det(do) = det(d 2 ). 



(37) 



Furthermore, the effective action (fT^|) was obtained by 
expanding the fermionic partition function Z[a^\ pertur- 
batively in the field a M , hence we get an additional pref- 
actor corresponding to the free fcrmion 



[0]-exp(i5 O ff[^]). 



(38) 



Hence, the bosonization formula with all prefactors in- 
cluded reads 

for mi [Aft, = 0] • N A N W J Da^DW^ 

5(dM dct(d 2 ) exp (iSbosonK, W»,A^ . (39) 

When setting A^ = everywhere, the functional integral 
over the bosonic fields clearly has to be equal to unity. 
This is also apparent from the interpretation of as 
a Lagrange multiplier and a M as a gauge degree of free- 
dom: integrating over W M constrains a M = d^A and the 
Fadeev-Popov ghost cancels any spurious thermodynam- 
ics that might come from this integration. Then, a M is 
constrained to be a purely longitudinal field and the pref- 
actor N A makes sure that integrating over it gives a finite 
result equal to unity. 

In 1 + 1 dimensions, the specific heat of the bosonic 
fields calculated in the manner of Equation (j3"4"|) without 
keeping track of the ghost contributions will match the 
fermionic specific heat, but this may be a coincidence: 
The contributions from the Fadeev-Popov ghost have to 
cancel the bosonic parts and the fermionic determinant 
is generated when expanding the effective fcrmion action 
in perturbation theory. 



VI. REFERMIONIZATION 

The usefulness of the bosonization procedure would 
be greatly increased if we could somehow express the 
original electron quantum field in terms of the collec- 
tive bosonic excitations. Such a refermionization pro- 
gram has been a great success in 1 + 1 dimensions [lj], 
as it allows us to study strong electron interactions be- 
yond perturbation theory. Two separate proposals for 
refermionization in 2 + 1 dimensions have been put forth 
[321 [33l | . We now show that these two approaches can 
be connected and we will demonstrate that a concrete 
expression for the fermion operator in terms of the field 
W,j, can indeed be constructed. 

In the context of topological insulators, the possibility 
of refermionization was discussed in Ref. (22|. We will 
very much follow this lead, though we have to make an 
important amendment concerning the physical interpre- 
tation of the bosonic fields first. 
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£(q) 




-> q 



FIG. 1: The shaded area represents the possible energies 
■E(q) = u|k + q| — u|k| for particle-hole excitations of massless 
Dirac fermions. The hole has momentum k. 



The first refermionization proposal by Luther [321 ] uses 
the so-called tomographic transformation [5ll If32j to map 
2 + 1 dimensional bosons with a linear dispersion rela- 
tion E(q) = u|q| to Dirac fermions. We will discuss this 
method in more detail in a moment, but our main con- 
cern right now is the following observation: It appears 
that we cannot identify these bosons with the charge or 
current density. Namely, if we look at the bosonized ac- 
tion for the fermions (fTT))) and use the identification (fTB|) 
for the current j 7 ' , we see that the correlation function for 
the current is essentially the Green's function ([2Tj]) of the 
square root (— D) 1 / 2 . But bosons with a linear disper- 
sion have an entirely different Green's function, namely 
that for the wave operator (— □). This discrepancy is 
also apparent if we consider the excitation spectrum of 
the Dirac fermions, i.e. the possible energies for creating 
an electron- hole pair. If an electron eigenstate with mo- 
mentum k is excited to a state with momentum k + q, 
the energy difference is 



E{q) = w|k + q| -«|k| 



(40) 



Unlike in 1 + 1 dimensions, the energy of this particle- 
hole excitation is not determined by the momentum q 
only; its possible values can be summarized in Figure [1] 
The extremal parts where the electron momentum k is 
parallel to q do have a linear dispersion, but other con- 
tributions to the particle-hole continuum disperse in a 
different fashion and do not correspond to bosonic fields 
with a linear dispersion relation. 

In other words, while Luther's construction yields an 
electron operator in terms of some bosonic fields, it ap- 
pears that these bosonic fields cannot be related to phys- 
ically relevant quantities in a straightforward manner be- 
cause their correlation functions differ significantly. As 
a consequence, this also means that the refermioniza- 
tion of the topological insulator in [22| suffers from the 
same problem of not allowing for a direct relation of the 
bosonic fields to the charge density. 

To make progress, we now consider the second 
refermionization proposal by Marino [33j . Essentially, 
this approach starts by considering the action (fTH)) . Af- 
ter integrating out the a^ field, one obtains an action 



expressed solely in the field W^. Using the calculation 
from Appendix DO in particular formula (L47|) . we obtain 



7T 2 (l + a 2 ) 



d 2 xdt 



ae 



^ X W„d,Wy 



(41) 



We have written the transversal component using W^ v = 
d^W v — d v Wn to highlight the connection to the Maxwell 
term and abbreviated the prcfactor of the Chcrn-Simons 
term with a = ||. The identification (|16p of the current 
j** with the curl of the field suggests that the creation 
of an electron corresponds to the creation of a vortex in 
the vector field W^. Hence .Marino proceeds to construct 
a vortex creation operator [33|, HH and identifies it with a 
fermion creation operator. However, these fermions need 
to be subjected to a transformation which is non-local 
in both time and space in order to obtain the standard 
Dirac fermions. 

We will not pursue this second path here. Rather, 
we wish to employ the tomographic transformation, but 
apply it to the physical field that encodes the charge 
density and current. The main idea taken from Marino's 
approach is that we have to apply a transformation that 
is non-local in both time and space. In particular, if we 
rescale the fields as 



W„ := (-D) 1/4 W M 



(42) 



the action becomes 



d 2 xdt 



1 



-as^W^-nf^dvWx] (43) 



up to an overall prefactor. It features a standard Maxwell 
term and some kind of "non-local Chern-Simons term" . 
Both terms have the same order of derivatives, namely 
quadratic. 

Now is a good time to apply the tomographic trans- 
formation. We will see that both the Maxwell and the 
Chern-Simons term, as well as their combination above 
transform very pleasantly. But first, let us review the 
mathematics of the tomographic representation, though 
our discussion will necessarily be brief. For additional 
details, we refer to Ref. [Fll |52|. 

The basic idea of the tomographic transform is to think 
of a field (f> in two-dimensional space not as a collection of 
amplitudes <^>(x) associated with points of space, but to 
think of them as amplitudes cj>(9, y) associated with lines 
of the two-dimensional space. A line is specified by an 
angle 6 £ [0,7r) indicating direction and a number y S M 
that indicates its signed distance to the origin, as illus- 
trated in Figure fVTl The amplitude associated with a line 
is essentially the sum of the amplitudes associated with 
the points on that line, i.e. we define a unitary transfor- 
mation 

(T<f>)(9, y) := / d 2 x S^y - n ■ x)^(x) (44) 
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term 



FIG. 2: A line (solid) is represented by the angle 9 of its 
normal vector and its distance from the origin y. The points 
x = (x\,X2) on this line are precisely those that fulfill the 
equation y — nx = 0. 



where n = h(9) = (cos9,s'm9) T is a unit vector perpen- 
dicular to the line and S 1 ^ 2 is a delta function 

/>oo 

5 1/2 (y-hx) := / dkk^e-^ie^-^+e-^-^) 
Jo 

(45) 

with the regulator r\ — > in the end. The choice of 1/2 for 
the exponent ensures that the transformation is unitary, 
so that for any two fields </>, \, we have 



d 2 x 0(x) X (x) = / dffdy (T(/>)(0,y)(Tx)(6,y) (46) 



where appropriate prefactors have been absorbed into the 
integration J dOdy over angles and distances. It was first 
noted by Radon [54| that this transformation can be in- 
verted, so that the field </>(x) can be reconstructed from 
its integrals over lines, though we will make no use of this 
fact here. 

The most striking feature of the tomographic trans- 
form is how it maps derivatives of the fields. In particu- 
lar, the Laplace operator becomes a simple square of the 
derivative d y independent of the angle 9 



T{V4>) = hd y {T<l>) 

t(a^ = dim. 



(47) 
(48) 



For vector fields like W = (Wi,W2) , it is convenient 
to define (spatially) longitudinal and transversal compo- 
nents as 



W L := h ■ TW = mTWx + n 2 TW 2 



(49) 



W T := h x TW = n 2 TW x - n 2 TW x (50) 
as this allows us to write (spatial) divergence and curl as 



T(V • W) = dyW 1 
T(V x W) = dyW 1 



(51) 
(52) 



Together with unitarity, these transformation rules im- 
ply the following remarkable expressions for the Maxwell 



rf 2 ., 



dOdy 



(dow 1 y - (dyW 1 y + (cw^ - d y w y\ (53) 

and the Chern- Simons term 

J d 2 xe^ x W tl d v W x = J d9dy 2W T (d W L -d y W ). 

(54) 

In particular, the combinations W T and Wql '■= 
(doW L — dyWo) are invariant under gauge transforma- 
tions Wfi — > W,x + d^A. To keep notation simple, we 
have used the same symbol W for both the original field 
and its tomographic transform here. 

The action (|4"3")) for the field can now be written in 
tomographic space as 

SfWjj = J dOdydt [(d Q W T ) 2 - (d y W T ) 2 + (W 0L ) 2 

-2aW T (d 2 y ~d 2 )^ 2 W QL \. (55) 

Shifting 

W T = W T + a(d 2 - O^-^Wol (56) 

absorbs the Chern-Simons term in the component Wql 
whose dynamics are trivial. We obtain 



S[W T ,W 0L ] 



d9dy 



W T (d 2 -d 2 )W T 

+ (1 - a 2 )W 0L ) 



(57) 



Discarding or integrating out the uninteresting Wql com- 
ponent, we see that the action for the W T field is pre- 
cisely the action needed for refermionization in 1 + 1 di- 
mensions. The equation of motion is simply the wave 
equation 



(d 2 d 2 )W T = 



(58) 



and we can decompose the field into left- and right- 
moving parts 

W T (6,y,t) = Wi e{t (9,y + t) + W lis u (8,y-t). (59) 

Hence, for every angle 9, we can apply the refermioniza- 
tion procedure familiar from 1 + 1-dimensional bosoniza- 
tion [l4| to obtain fermion fields as a normal ordered 
version of 

-ip±(0,y) = Cexp (iViTW leit/light (8,y,t)) 6 e (60) 

where C is some constant and Og refers to Klein factors 
that are needed to obtain the correct anticommutation 
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relations. Undoing the tomographic transform yields a 
massless field of 2 + 1 dimensional Dirac fcrmions again, 
though we refer to Ref. [HI, [f52j for the details on this 
construction. 

To summarize, the tomographic representation allows 
us construct a fermion operator ip from the field with 
the same Green function as the 2 + 1-dimensional mass- 
less Dirac fermion. From this, we can see that the cur- 
rent tpj^ip of this new fermion has the same correlation 
functions as the current j M = e^ vX d u W\ of the original 
fermion field. Hence, we propose to identify the original 
fermion field with the newly constructed fermion. 



effective bosonic action (fl9|) is a strong indication that 
the bosonic surface theory actually describes fcrmionic 
excitations. The Klein factors in the refcrmionization 
proposal (|60p are yet to be constructed explicitly in terms 
of the bosonic field W M , but they may be related to this 
statistical transmutation. 
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for collaborating with us in the initial stage of this 
project, Daniel Scherer and Mirco Milletari for help- 
ful discussions, Norman Metzner for feedback on the 
manuscript and BMBF for financial support. 



VII. SUMMARY AND DISCUSSION 



Appendix A: Field Integral Calculations with the 
Chern-Simons Term 



Guided by the principle of anomaly cancellation be- 
tween bulk and surface, we have constructed a topologi- 
cal field theory for the TRS topological insulator in three 
space dimensions. We have incorporated the electron 
dynamics into the bosonic surface theory by consider- 
ing the effective electromagnetic response obtained from 
the fermionic Dirac Hamiltonian. This bosonic theory 
does not reproduce the fermionic thermodynamics quan- 
titatively, but we have shown that this discrepancy is 
accounted for by normalization factors of auxiliary func- 
tional integrals. Using the tomographic transform, we 
have proposed a representation of the fermionic field in 
terms of the bosonic field, though the non-locality of the 
bosonic action forces this representation to be highly non- 
local as well. 

Our approach is very much inspired by the construc- 
tion of a topological field theory for the TRS topological 
insulator in three dimensions in Ref. [22:]. In contrast to 
Ref. [22[, however, we have argued that the TRS insulator 
is described, perhaps counterintuitivcly, by a surface the- 
ory that explicitly breaks time-reversal symmetry, even 
though the combination of bulk and surface is again time- 
reversal invariant. Also, we have found that the effective 
dynamics of the charge density on the surface are highly 
non-local in both space and time, as can be seen by ex- 
amining the spectrum of particle-hole excitations. 

Similar to the construction of topological field theories 
in Ref. [2(J, the method of functional bosonization has 
allowed us to derive systematically a field theory from 
the response to an external field. However, our main 
guiding principle is not just the electromagnetic response, 
but the cancellation of quantum anomalies, which has 
allowed us to predict the existence of surface states from 
the topological bulk theory. Moreover, we were able to 
identify the appropriate degrees of freedom needed for a 
topological field theory of the quantum spin Hall effect 
by discussing the axial anomaly in Appendix IBI 

We believe that the parity anomaly may have addi- 
tional significance in the context of refermionization. The 
Chern-Simons term is known to induce a statistical trans- 
mutation [55l . [56j , mapping bosons to fermions and vice- 
versa. Hence, the presence of a Chern-Simons term in the 



We now introduce some tools which facilitate the 
calculation of functional integrals involving the 2 + 1- 
dimensional Chern-Simons term. The key convenience 
is that there exist [43j two orthogonal projection matri- 
ces P + , P~ that project a vector onto its (spacctimc) 
transversal components 



(P+ +P-) llu a v = S„ u - 



(Al) 



but are also related to the Chern-Simons [6l| term via 

a fi \q\(iP + - iP~)^a v = e^a^a^. (A2) 

These two matrices are readily obtained from the "Chern- 
Simons matrix" 



• Ha 

m 

as the relations 

* f * f * 



r3 r~< { ■ f > { ' f ■ 



allow us to write 
1 



P, 



(c'±c^), {p ± y^p ± , p+p 



(A3) 

(A4) 
(A5) 



= 0. 

(A6) 



With these tools, we can write the surface action (fT 
in Euclidean momentum space as 

S[o M , W„] = [ ~ ^(-q)\q\( t P+ - iP-)^W v {q) 



+ _o„(-g)|5|((l + ia)P+ + (1 - ia)P-)^ v a„(q) 



(A7) 



using the notation a 
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For instance, it is now straightforward to integrate out 
the field a M and obtain the action expressed in terms of 
the field alone. Namely, the linear shift 

+ Z (t^~ P+ - T^— p ~) W » ( A8 ) 



7r V 1 + ia 1 — ia 



fJjV 



will complete the square and leave the residual action 

TT i (l + or) ' 
Wu\q\ {P+ + P~ - ia(P+ - P-)) ^Wu (A9) 
in Euclidean momentum space. 



Appendix B: Axial Anomaly and the Quantum Spin 
Hall Insulator 



We now want to come back to the discussion of quan- 
tum anomalies from Section |TT] and explain the edge 
states of the TRS topological insulator in two space di- 
mensions in terms of the axial anomaly [35l |41| . For 
concreteness, we will focus on the special case of the 
quantum spin Hall insulator [57j], and assume that the 
electrons in the planar bulk carry a conserved quantum 
number, namely the spin in perpendicular direction, S z . 
This model system is representative for TRS topologi- 
cal insulators (58|, but we admit that it is not entirely 
clear how to find a substitute for the spin current in the 
general case. Still, we wish to argue that it beautifully 
captures the topological essence of the situation. 

The quantum spin Hall insulator can be understood 
as two opposing copies of the quantum Hall effect: The 
spin-up electrons experience a Hall effect with conduc- 
tivity <J xy , while the spin-down electrons are subject to 
a Hall conductivity —o~ X y, i.e. they move in the opposite 
direction. Time reversal corresponds to interchanging 
the spin up and spin down directions, and the system is 
therefore invariant under time-reversal. Since the elec- 
trons move in opposite directions, the total electric Hall 
current vanishes, but we can observe a spin Hall current. 

Again, the key feature of the quantum spin Hall insu- 
lator is that the edge of a finite sample features gapless 
electron states. Since spin is conserved, we could use 
our knowledge of the quantum Hall effect and the chiral 
anomaly to predict that the edge has two counterpropa- 
gating modes: the spin-up electrons moving in one and 
the spin-down electrons moving in the other direction 
along the edge. However, anticipating a generic topo- 
logical insulator where spin is not conserved, we want 
to answer the following question: How can we obtain 
the combined edge state from a quantum anomaly with- 
out separating the electron into independent spin com- 
ponents? 

The main observation is that the combination of the 
two edge modes is equivalent to the 1 + 1 dimensional 
Dirac fcrmion, which is known to suffer from the axial 



anomaly [35|,|41|]. Essentially, the axial anomaly is a com- 
petition between the charge current = ipj^x/j and the 
spin (or axial) current = tpj^j^tp: while the classical 
action for the Dirac fermion suggests that both currents 
satisfy the continuity equation, the quantum theory ac- 
tually predicts that the currents cannot be conserved si- 
multaneously, we can choose between either dg (j^) = 
or ^hUb) = 0; but we can never have both 62|. But 
this necessity to choose also makes it clear that we can- 
not predict the edge modes from considering the charge 
current or the spin current alone, because each of them 
can be conserved individually and we would not observe 
a defect. 

From this point of view, we must consider an effec- 
tive bulk action S^A^, B^] that depends on both the 
electromagnetic gauge field and a (fictitious) axial 
gauge field B^ so that we can probe both the charge cur- 

}M _ SS[A,B] 
1 ~ <5A M 
SS[A,B] 
SB U 



rent via j 7 



Jb 



and the axial (spin) current via 

Essentially, we think of it as a generat- 
ing functional for the two currents. Invariance of the 
action under gauge transformations A^ — s- A^ + d^A or 
B^ — » B^ + <9 M means that the charge current, respec- 
tively the axial current is conserved. The effective action 
is easily engineered from the Chern-Simons theory of the 
quantum Hall effect |34| . Representing the spin- up elec- 
trons by a collective field a and the spin-down electrons 
by a collective field b, we couple the electromagnetic field 
with the same sign, but the axial field with opposite signs 



S[a,M M ,B M ] = J- / dI xdt 



1 



1 



2tt 

^VcUa + e^A^ax + s^ x B^a x 
e^%dM + e^A^bx - e^B^dM 



(Bl) 



Integrating out the matter fields, or simply noting that 
the spin-up electrons basically see a total gauge potential 
of A^ + B^ whereas the spin-down electrons see a total 
gauge potential of — A^ + B^, we obtain the effective 
action 

S cS [A^B^ = a xy [ d 2 xdt Is^A^Bx+e^B^A 



(B2) 

As in our previous discussions, for periodic boundary 
conditions, the effective action is invariant under gauge 
transformations of the fields A^ and B^. In fact, integrat- 
ing partially, we can write the action in many different 
ways 

Seff[A„, £y = 2a xy J d 2 xdt e^BpdvAx (B3) 



= 2a xy J d 2 xdt e^ApdvBy 



(B4) 



so that the action is manifestly gauge invariant in one 
field, A^ in the first equation, because the action only 
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depends on its curl, and gauge invariant in the other field, 
in the first equation, because of the periodic boundary 
conditions. In the second equation, the manifestation of 
the gauge invariance is reversed. 

Now, if we try to restrict the bulk action to a bounded 
domain, say to the region x\ > 0, we can choose between 
many different forms like $B% . $B3\i and ([52]) for the bulk 
action. While puzzling at first, this choice is precisely the 
ambiguity that the axial anomaly allows: we may choose 
to conserve either the charge or the spin current! When 
choosing the form (|i33j) where gauge invariance of the 
electromagnetic field is manifest, the charge current will 
be conserved, but the spin current will be anomalous and 
must be cancelled by a corresponding edge mode 

d,lU = (B5) 
«,buik = 2a xy S(x 1 )s 1 ^d,A x . (B6) 

Similarly, the form (|i34[) must have an edge mode that 
cancels an anomalous charge current while the spin cur- 
rent is conserved. Mixed forms are also possible. 

To summarize, we have shown that the bulk action(s) 
(|_B2I) predicts the edge states of the quantum spin Hall 



insulator via the axial anomaly. In particular, it correctly 
reproduces the ambiguity inherent in the axial anomaly 
and shows that one must consider both charge and spin 
current in combination to predict the edge state. Again, 
it is straightforward to construct a corresponding topo- 
logical field theory via bosonization as presented in Sec- 
tion nu 

With regards to the ambiguity of the axial anomaly 
and the physical meaning of the axial current, in a con- 
densed matter physics setting, it seems likely that the 
charge current is conserved at the edge and that the 
spin current must carry the defect. However, such a 
conclusion cannot be drawn from topological considera- 
tions alone. The origin of the choice inherent in the axial 
anomaly is the need to regularize the quantum theory at 
high momenta, and the availability of several choices for 
doing so. In a condensed matter setting, a lattice regu- 
larization is part of the theory, and the details depend 
on the physical system at hand. Likewise, whether the 
axial current corresponds to a spin current or to some 
other quantity related to time-reversal symmetry should 
depend on the details of the system as well. 
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